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In this paper we study the sound speed - c 2 a , which is directly related to the classical perturbations 
- of the dynamical dark energy (DE) , especially with an equation of state crossing the cosmological 
constant boundary in details and show its implications on Cosmic Microwave Background (CMB) 
Anisotropy. With the present observational data of CMB, Type la Supernova (SNIa) and galaxy 
clustering, we perform a global analysis to constrain the sound speed of dark energy, using the 
Markov Chain Monte Carlo method. We find that the sound speed of dark energy is weakly con- 
strained by current observations thus the futuristic precision measurements of CMB on a very large 
angular scale (low multipoles) are necessary. 
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I. INTRODUCTION 

The analysis of the redshift-distance relation of SNIa revealed that our Universe is currently accelerating [l], E| 
which has been confirmed by the current h igh quality observations of the CMB [3(, the Large Scale Structure (LSS) of 
galaxies [S El, El 0] i an d the SNIa d H [H EH • One possible explanation for this phenomenon is that this acceleration 
is attributed to a new form of energy, dubbed Dark Energy, which recently dominate the energy density of the Universe 
with negative pressure and almost not clustering. The nature of DE is among the biggest problems in modern physics 
and has been studied extensively. A cosmological constant, the simplest DE candidate whose equation of state w 
remains — 1, suffers from the well-known fine-tuning and coincidence problems. Alternatively, dynamical dark energy 
models with rollin g^ sc alar fields have been proposed, such as Quintessence 12 , H H3] , Phantom [HI , Quintom fl6j 

and K-essence IS SEE!. 

Given our ignorance of the nature of dark energy, the cosmological observations play a 



crucial role in our understan ding of DE. There are many studies on DE both theoretically and phe nomenologically in 

the literature [U El El El ES El, El El El, El EI M, El El ES ES ES El El EI EI EH, Si, El E3, EE EI El 
" HE!- 

If we regard the mysterious DE as a cosmic fluid, there will be perturbations when w ^ — 1. In the framework of 
the linear perturbations theory, the DE perturbations can be fully characterized by its equation of state, w = p/p, 
and its sound speed, c 2 s = Sp/Sp. For the DE per turbations it is well-defined in the region of Quintessence, w > — 1, 



and Phantom, w < — 1. In our previous works H 23[ based on Quintom dark energy model, we have proposed a new 
method to handle the DE perturbations in the region closed to w — — 1, assuming the sound speed of Dark Energy 
is equal to unity, (? s = 1. We found that the constraints on the DE parameters are relaxed and the parameter space 
are enlarged dramatically when the DE perturbations are included in the analysis. The effects of perturbations on 
studying the feature of DE are of great significance. 

Since the sound speed of DE merely affects the evolution of perturbation, it has not been investigated extensively 
in the literature. Recently the observational implications of the sound speed on CMB and LSS of some DE models 
have been discussed: for exam ple, for K-essence 51, EH, condensation of Dark Matter [53] and the Chaplygin gas 
H El, El, El HI • In Refs. [H El EH the authors tried to use the current observational data to constrain the sound 
speed of DE. However, they just considered the equation of state of DE to be constant. In this paper we consider 
the dynamical DE model with time varying w and extend our previous works on the sound speed of DE models and 
allow an arbitrary constant sound speed in our calculation. Combining the present observational data, such as CMB, 
LSS and SNIa, we discuss the possible constraints on the sound speed of DE. Our paper is organized as follows: In 
Section II we analyze the feature of the sound speed of dark energy models paying attention to the case when the 
equation of state gets across w = — 1; In Section III we describe the method and the datasets we used; In Section 
IV we present our results derived from a global analysis using Markov Chain Monte Carlo (MCMC) method; Finally, 
Section V contains a discussion of the results. 



II. FEATURE OF SOUND SPEED 



A. Single Perfect Fluid 



Working in the conformal Newtonian gauge, one can easily describe the DE perturbations as follows [62J: 

6' = -(l + w)(0-3$') -m{cl-Wi)6 , (1) 
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9' = -ft(l 3w)9 —6 + k 2 ^ + *) , (2) 

where the prime denotes the derivative with respect to conformal time, c 2 = Sp/8p denotes the sound speed of DE 
models, 5 and 9 are the density and velocity perturbations of DE models respectively, w = p/p is the equation of 
state, and 4> and \I> represent the metric perturbations. 

Let us discuss the behavior of the perturbations of single perfect fluid when w crosses the boundary w = — 1 in the 
conformal Newtonian gauge, following Firstly, we consider the DE as a single barotropic fluid. The adiabatic 
speed of sound, c 2 a , is purely determined by the equation of state w, 

Sp = S p = ctSp =(w- m ™ +w) ) Sp . (3) 

Whenever the equation of state closes to the boundary w = — 1 and w'\ w —-i ^ 0, the adiabatic speed of sound, c 2 , 
will be divergent due to the existence of the term w'/(l + w). Apparently, the perturbations of this system become 
unstable. 

However, in non-barotropic fluids, for example in most scalar field models, this simple relation in Eq.Q between 
the equation of state and the sound speed breaks down, because of the intrinsic entropy perturbations, and we have 
the more general relation 

In this case, if c 2 ^ c 2 , the intrinsic entropy perturbation T will be induced 6(| 63 1: 

P \P P J 

Whereas the adiabatic sound speed, c 2 , and T are scale independent and gauge invariant quantities, while c 2 can be 
neither. In order to construct the gauge invariant c 2 we use a helpful transformation [63j: 

8 = 5 + 3H(l + w)-z . (6) 

This transformation relates the gauge invariant, rest frame density perturbation, 8, to the density and velocity 
perturbations in a general frame, S and 9. Using Eqs. (|5l6p . we can rewrite the pressure perturbation in a general 
frame, 8p, in terms of the rest frame sound speed, c 2 , 

5p = c 2 s 5p + 3W(1 + w)(c 2 s - c 2 a ) P -^ . (7) 

For example single scalar field models with standard kinetic terms always have c 2 = 1. Moreover, the Eqs. (|l|2p can 
be rewritten as: 

Q 

5' = -(1 + w)(9 - 3$') - m(c 2 s - w)S - 3H(w' + 3ft (1 + w)(c 2 s - w))-^ , (8) 

k 

9' = -H{l--Zc 2 )9 + k 2 {^- + m) . (9) 

1 + w 

In this case if we require the DE equation of state can cross the boundary w = — 1, the DE perturbations are still 
divergent (see appendix for details). 



B. Two-field Quintom Dark Energy Models 



As discussed above, the DE models, whose equation of state can cross the boundary w 
of freedom to keep the whole system stable 



64|. 



with two components, one is Quintessence-like, 
analyze the perturbations of this DE system. 



1, need more degrees 
As an example, we consider one type of Quintom DE models 

"TSUI HI and 



T < w < 1, and the other is Phantom- like, w < — 1 
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The most notable advantage, in this case, is that it is unnecessary to let the total equation of state, lutot, of 
this Quintom system cross the boundary w = — 1 by forcing the equation of state of each component to cross the 
boundary. For the Quintessence- like component, — 1 < w < 1, while for the Phantom- like component, w < — 1, 
during the evolution, without any crossing behavior. For each component of the Quintom model the density and 
velocity perturbations, Si and 8i, still satisfy the Eqs. (|8l9p . Because the equation of state for each component needs 
not cross the boundary, the perturbations of each component will be stable and all the variables will be well defined. 
Furthermore, all these corresponding variables of the whole Quintom DE model can be combined with each component: 

(10) 
(11) 
(12) 



Ptot 


= Pi + P2 , 






Ptot 


= Pl+P2 , 

W\pl + W 2 p2 






Wtot 


Pi + p2 






<5tot 


pih + P2S2 
P1+P2 






Sptot 


= Spx + Sp 2 , 






Otot 


(l + tui)pi(9i - 


Ki- 


- w 2 ) P2O2 


(1 + wi)pi - 


Ki- 


V W 2 ) P2 



(13) 
(14) 

(15) 

Using these Eqs. (Til)l -[l"4"l). these variables of Quintom DE models are well defined except the velocity perturbations 

0tat- 

As long as w to t gets across the boundary, w to t — > — 1, the denominator of Eqs. (fl"5|) goes to zero, [(1 + wi) p~\ + 
(1 + W2) P2] — * 0. Consequently, the total velocity perturbations, 6 to t, become divergent unless these two independent 
perturbations are equal all the time, #tot — 8\ =62- During the evolution, it is impossible to keep the velocity 
perturbations of each component equal, 9\ =02, ah the time. However, the physically meaningful velocity perturbation 
which is relevant to the CMB observations is V = (1 + w)8, not 9 alone 22, 66|. Eq. ([T5|) can now be written as: 



,. V1P1+V2P2 ( , R , 

Hot = — _ , _ • (16) 

Pi + P2 

Using this Eq.(fT6|). the divergences disappear when the total equation of state crosses the boundary w = — 1. And 
then all the variables are well defined and the perturbations of Quintom dark energy system are stable in the whole 
parameter space of the total equation of state, adding more degrees of freedom. 

Let us move to the discussion of the sound speed of the whole Quintom dark energy system. Physically we can 
use the independent sound speed of each component to describe the whole system. However, the present constraints 
on the sound speed of DE are so weak that considering to study the two independent sound speed is not justified at 
present. For simplicity we investigate the effective sound speed, c? s eS [6^ |: 

.2 = (1 + w i) gl + g s,2 (1 + w 2 ) P2 

Cs ' cff ~ (l + u>i)pi + (l + wa)fr ' 1 ' 

It is apparent that this effective sound speed suffers from the same problem with the velocity perturbations # to t- In 
order to settle this divergence of the effective sound speed we fix the sound speed of each component to be equal, 
c 2 s 1 = (? s 2 = c? s . The effective sound speed of the whole Quintom system should be a constant during the evolution: 

S 2 _ C 2 S (1 +W 1 )p 1 +C 2 S (1 + W 2 )p2 _ , 2 _ -2 _ -2 ns s 

C s — C s.l — C s,2 ■ 



^s.eff — T\ 1 \ ^ 1 n 1 „ ,, \ — — — °s,l — ^s,2 



(1 +W\)p\ + (1 + w 2 )p 2 

The total equation of state of Quintom system has nothing to do with the effective sound speed now 



In our previous works [22j,|23| we considered the Quintom dark energy model with the effective sound speed c 2 scS = 1. 
In this paper we relax this limitation and assume an effective arbitrary constant sound speed of Quintom DE model, 
c 2 eS = c 2 . The effects of this effective sound speed on CMB power spectrum will be studied in the next subsection. 



C. Implications on CMB Anisotropy 



In Fig|T| we illustrate how the CMB temperature anisotropies of four different DE models charactrized by different 
equations of state w change on large scales, for different constant sound speed c 2 . Recently we know that our Universe 
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is dominated by the dark energy component. This means the sound speed of DE can only affect the CMB power 
spectrum at the very large observable scale via the late Integrated-Sachs- Wolfe (ISW) effect which describes the 
perturbations induced by the passage of CMB photons through the time evolving gravitational potential wells during 
DE domination 5^, 6(J 6l|, [67 1. The explicit effect on CMB power spectra caused by the sound speed of DE appears 
at the large angular scale / < 20. However, as the behavior of the equation of state is close to the Cosmological 
Constant boundary w = — 1, these differences get smaller and even disappear. The different effects of the two DE 
models in the bottom panels are much clearer than the ones in the top panels which are closer to the boundary. 
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FIG. 1: The effects on CMB temperature power spectra of four different Dark Energy models charactrized by different equation 
of states w, for different constant Dark Energy sound speed c? a . On the left two panels, the bottom Black solid line is for the 
sound speed of (? s = f .0 and the top Magenta dash dot dot line is with c? 3 — 0. In between the sound speed is increasing from 
top to down with c? s = 0.01,0.05,0.2. On the right two panels, the top Black solid line is for the sound speed of c? a = 1.0 
and the bottom Magenta dash dot dot line is with (? a = 0. In between the sound speed is decreasing from top to down with 
c 2 3 = 0.2, 0.05, 0.01. 



III. METHOD AND DATASETS 

During our calculation we choose the commonly used parametrization of the DE equation of state as p38j : 

wde(o) = wo + wi(l - a) , (19) 

where a = 1/(1 + z) is the scale factor and w\ — —dw/da characterizes the "running" of the equation of state. We 
modify the current publicly available codes like CMBFAST \fX^ and CAMB [70| to allow the DE equation of state 
to cross the boundary w = —1. For the parametrization of the equation of state which gets across -1, we introduce 
a small positive constant e to divide the full range of the allowed value of w into three parts: 1) w > — 1 + e; 2) 
— 1 + e > w > —1 — e; and 3) w < —1 — e. We neglect the entropy perturbation contributions, for the regions 1) and 
3) the equation of state does not get across —1 and perturbations are well defined by solving Eqs. (|ll2[ ). For the case 
2), the density perturbation 5 and velocity perturbation 9, and the derivatives of S and 9 are finite and continuous 
for the realistic Quintom Dark Energy models. However, for the perturbations of the parameterized Quintom there 
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is clearly a divergence. In our study for such a regime, we match the perturbation in region 2) to the regions 1) and 
3) at the boundary and set: 



5' = , 0' = 



(20) 



In our numerical calculations we have limited the range to be |e| < 10 5 and we find our method is a very good 
approximation to the multi-field Quintom DE model. For more details of this method we refer the readers to our 



previous companion papers 



In this study we have implemented the publicly available Markov Chain Monte Carlo package CosmoMC [71| , which 
has been modified to allow for the inclusion of DE perturbation with the equation of state getting across —1. We 
assume purely adiabatic initial conditions and a flat Universe. Our most general parameter space is: 



(u>6, u c , Q s , t, c 2 , w , wi, n s , log[10 10 A s ]) 



(21) 



where lo^ = Vt^h 2 and lo c = ft c h 2 are the physical baryon and Cold Dark Matter densities relative to the critical 
density, S is the ratio (multiplied by 100) of the sound horizon to the angular diameter distance at decoupling, 
t is the optical depth to re-ionization, c 2 is the effective sound speed of Dark Energy, A s and n s characterize the 
primordial scalar power spectrum. For the pivot of the primordial spectrum we set k s o = 0.05Mpc _1 . Furthermore, 
we make use of the Hubble Space Telescope (HST) measurement of the Hubble parameter Hq = lOOh km s _1 Mpc -1 
[72| by multiplying the likelihood by a Gaussian likelihood function centered around h = 0.72 and with a standard 
deviation a — 0.08. We also impose a weak Gaussian prior on the baryon density Vlbh 2 = 0.022 ± 0.002 (1 a) from 
Big Bang Nucleosynthesis [73j . Simultaneously we will also use a cosmic age tophat prior as 10 Gyr < to < 20 Gyr. 

In our calculations we have taken the total likelihood to be the products of the separate likelihoods of CMB, LSS 
and SNIa. In the computation of CMB we have included the three-year WMAP (WMAP3) data with the routine 
for computing the likelihood supplied by the WMAP team For LSS information, we have used the 3D power 
spectrum of galaxies from the SDSS [B| and 2dFGRS jij. To be conservative but more robust, in the fittings to the 
3D power spectrum of galaxies from the SDSS, we have used the first 14 bins only, 0.015 < fc c ff < 0.1, which are 
supposed to be well within the linear regime [6( ■ In the calculation of the likelihood from SNIa we have marginalized 
over the nuisance parameter [3]. The supernova data we used are the "gold" set of 157 SNIa published by Riess et 
al in Ref.fi]. 

For each regular calculation, we run 8 independent chains comprising of 150, 000 — 300, 000 chain elements and 
spend thousands of CPU hours to calculate on a supercomputer. The average acceptance rate is about 40%. We test 
the convergence of the chains by Gelman and Rubin criteria 75[ and find that R — 1 is of order 0.01 which is more 
conservative than the recommended value R — 1 < 0.1. 




0.5 -10 12 



FIG. 2: Constraints in the (u;o, <? a ) and (wi, c? s ) planes at 68% (dark) and 95% (light) C.L. from a combined analysis of CMB, 
LSS and SNIa observational data together. 
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IV. GLOBAL FITTING RESULTS 

In this section we mainly present our global fitting results of the dark energy parameters wq, w\ and the effective 
sound speed of dark energy cf. Firstly, in Fig [2] we show the constraints on the effective sound speed c 2 from a 
combined analysis of the CMB, LSS and SNIa observational data. The likelihood contours of (wq, c 2 ) and (w\, c 2 ) 
at 68% and 95% C.L. are nearly vertical lines. The dark energy parameters, as well as other cosmological parameters 
which we do not list here, are almost independent of the effective sound speed of DE. There is nearly no constraint on 
the effective sound speed from the present astronomical data, namely, the current observations are still not sensitive 
to c 2 . Thus the futuristic precision measurements of CMB on a very large angular scale (low multipoles) are necessary. 



1.5 



0.5 
-0.5 



-1.5 

-2 

-1.4 -1.2 -1 -0.8 -0.6 

w o 

FIG. 3: Two-dimensional constraints on the DE parameters wo and wi from a combined analysis of CMB, LSS and SNIa data 
at 68% and 95% C.L. with c? s = 1 (Black Dashed lines) and varying (Blue Solid lines). The Magenta Dotted lines stand for 
wq — — 1 and Wo + w i = — 1. The Roman letters denote different dark energy models. See text for details. 

On the other hand, we consider the effect of the effective sound speed on the constraints of DE parameters. In Fig[3] 
we compare the constraints on the dark energy parameters wq , wi with c 2 = 1 (Black Dashed lines) and with arbitrary 
c 2 (Blue Solid lines). As expected, the constraints on wq, Wi almost unchange. For the model with c 2 = 1, we find 
that w = -1.05to.\ltoMi and w i = 0.53±g;||±?:^|, meanwhile, w = — 1.06±g;i|±o;lI and w i = °- 54 ±o:S±i'.52 for thc 
c 2 varying models. We cut the parameter space of wo — w\ plane into six parts by the line of Wo = — 1, wq + Wi = —I 
and W\ — 0. Part III is for quintessence-like models, namely, the equation of state remains greater than —1 regardless 
of cosmic time, say, w > — 1 for past, present and future. Correspondingly, part VI is for Phantom-like models. Part 
I,II,V and IV are all for Quintom-like models. For the models lie within part I and IV, their equations of state have 
crossed over -1 till now while the EoS of the DE models in part II and V will cross -1 in future. 

V. SUMMARY 

In this paper we have studied the features of the sound speed of DE in detail and have used the present observational 
data to constrain the effective sound speed of dark energy c 2 and the equations of state Wo, w\. 

If the sound speed is smaller than zero the system is unstable, due to the divergent classical perturbations. If we 
regard the DE models as single barotropic fluid or non-barotropic fluid, when the equation of state w is close to the 
boundary w = — I and w'\ w —-\ ^ 0, the classical perturbation of the DE system will be divergent. In order to settle 
this problem it is necessary to add more degrees of freedom. As an example, we analyze the stability of Quintom DE 
models with two components. We assume that the two components have the same sound speed to avoid the intuitive 
divergence. 

Using the Markov Chain Monte Carlo method, we preform a global analysis of the c 2 and Wo, W±. We find that 
the current astronomical data have nearly nothing to do with constraining the effective sound speed of Dark Energy 
system. The constraint on the sound speed of DE is very weak. The futuristic precision measurements of CMB on a 
very large angular scale (low multipoles) are necessary. 
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APPENDIX: THE PROOF OF THE NO-GO THEOREM 



In this appendix we present the detailed proof of the "No-Go" Theorem which forbids the equation of state parameter 
of a single perfect fluid or a single scalar field to cross the —1 boundary. 

Theorem: For theory of dark energy (DE) in the Friedmann-Robertson- Walker (FRW) universe described by a 
single perfect fluid or a single scalar field cj> with a lagrangian of C — C^jd^d^fi), which minimally couples to 
Einstein Gravity, its equation of state w cannot cross over the cosmological constant boundary. 

Proof: Let us consider the case of fluid firstly. Generally, a perfect fluid, without viscosity and cannot conduct 
heat, can be described by parameters such as pressure p, density p and entropy S, satisfying the equation of state 
P = p(p> According to the properties of fluid, single perfect fluid can be classified into two kinds of form, dubbed 
as barotropic and non-barotropic. 

If the fluid is barotropic, the iso-pressure surface is identical with the iso-density surface, thus the pressure only 
depends on its density in form of p = p(p). From Eq.([3]), we can see that the sound speed of a single perfect fluid is 
apparently divergent when w crosses —1, which leads to instability in Dark Energy perturbation. 

If the fluid is non-barotropic, the pressure generally depends both on its density and entropy, p = p(p, S) . The 
simple form of the sound speed defined in Eq.([3]) is not well-defined. From Eq.lfJJ while taking gravitational gauge 
invariance into consideration, we can obtain a more general relationship between the pressure and the energy density 
as follows, 

5p = c 2 s 5p , (22) 
where the hat denotes the gauge invariance. Eq.([6]) can be written in a more explicit form as: 

6p = 6p + 3H{p + p)-jL , (23) 

and correspondingly, gauge-invariant perturbation of pressure turns out to be 

5p = 5p + 3Hcl(p+p)^ , (24) 

where p and p are background energy density and pressure respectively, while 9 is the perturbation of velocity as 
defined in the main text. The gauge-invariant intrinsic entropy perturbation T can be described as, 

T = —_(5p- c 2 Jp) = —_(6p- c 2 Jp) . (25) 
wp wp 

Combining Eqs. (|22H24|) . we obtain the following expression, 



Sp = c^dp 



,. , , mpo{i + w){ci-<? a ) 



k 2 



= ^ + ^^±^.^±^ m+ s^. (26) 

From the Dark Energy perturbation equation Eq.Q, one can see that 9 will be divergent when w crosses —1, unless 
9 satisfies the condition 

e w > = k 2 % . (27) 



So we check what will occur when the condition is satisfied. 
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By substituting the definition of adiabatic sound speed c 2 and the condition Eq. (|27|) into Eq. ([24|) , we obtain Sp = 0. 
Note that T = ^=(5p — c 2 Sp), it is obvious that due to the divergence of c 2 at the crossing point, we have to require 
Sp = to maintain a finite T. So we come to the last possibility, that is, dp = and dp — 0. From Eqs. (|23l24p . this 
case requires that Sp = —c 2 3H P^ i 1 +" J ) anc ! Sp = — ■ in P ^-+ w ) ; anc j ^hus Sp = c 2 a Sp. It returns to the case of adiabatic 
perturbation, which is divergent as mentioned above. 

Finally, from the analysis of classical stability, we demonstrate that there is no possibility for a single perfect fluid 
to realize w crossing — 1. For other proofs, see 64. |65||. 

In the following part we will discuss the case of a single scalar field. The analysis is an extension of the discussion 
in 22j . The action of the field is given by 

3= [ d^x^C^d^d^) , (28) 



where g is the determinant of the metric g^. To study the equation of state (EOS) of DE, we firstly write down its 
energy-momentum tensor. By definition that S g ^ u S = — J d 4 x^^-T p ' v Sg^ u , one can get 

Sg^v oX 

where X = ^d^d^cf). In FRW universe, we have the metric in form of g^ — diag(l, —a 2 , —a 2 , —a 2 ). In the framework 
of Einstein Gravity, we can neglect the spatial derivatives of cj> and rewrite X as X = \4> 2 , where the overdot denotes 
the derivative with respect to cosmic time t. So T^ 11 has the matrix form of diag(p, \), comparing with the 

fluid definition: T" v = {p+p)U^U u - pg» v (or T£ = (p + p)U"U v - pS£) where V is the four-velocity of the fluid. 
Thus one can get: 



p = -n = c , (30) 

p = n = -sic + ^d°<pd <j> , 

= 2Xp tX ~P, (31) 



where " t x" stands for "^f". 

Using the formulae above, the equation of state w is given by: 

w = — = = — 1 H . (32) 

p 2Xp. x - p 2Xp iX - P 

This means that, at the crossing point t* , Xp^x\t* = must vanish. Since w needs to cross —1, it is required that 
Xp t x changes sign before and after the crossing point. That is, in the neighborhood of t* , (t* — e, t* + e), we have 

Xp >X \t*-e-X P ,x\t*+e<0 ■ (33) 

Since X = ^<j> 2 is non-negative, the Eq. (|3"3"|) can be simplified as p.x\t*-e • P,x\t*+e < 0. Due to the continuity of 
perturbation during the crossing epoch, we obtain p,x\t* — 0. 

Here we consider the perturbation of the field. We calculate the perturbation equation with respect to conformal 
time rj in form of: 



c 2 y 2 u-[— + 3c 2 s (H' -H 2 )}u = , (34) 



where we define 



S'2 



P,x\ , (35) 



and prime denotes and "H — -4r-". 

When expanding the perturbation function u by Fourier transformation, one can easily obtain the dispersion 
relation: 

c 2 = c 2 fc 2 -^-3c 2 (H'-H 2 ) , (36) 
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with c 2 defined as To make the system stable, we need c 2 > 0. Note that at the crossing point we require that 
p_x = and p t x will change its sign during crossing, one can always find a small region where c 2 < unless p t x also 
becomes zero at the crossing point with the similar behavior of p,x- Therefore, the parameter z will vanish when 
crossing. 

At the crossing point, TC' — Tt 2 in the last term in the expression of w 2 is finite, and when assuming that the universe 
is fulfilled by that scalar field, it turns out to be zero. 

Since at the crossing point, we have p t x = 0, z — 0. Consequently, if z" ^ at the crossing point the term will 

be divergent. Note that, even if z" = this conclusion is still valid. In that case, 2— = due to the L'Hospital 
theorem. Since z is a non-negative parameter, z = is its minimum and at that point z' must vanish, 2-r- is either 

divergent or equal to ^77- where z" is also equal to zero as we discussed before. Along this way, if we assume the 
first (n — l)-th derivative of z with respect to r\ vanishes at the crossing point and z^ ^= 0, which can be applied 

to an arbitrary positive integer n, we can always use the L'Hospital theorem until we find that ^- = jt-2) 1 which 
will still be divergent. Therefore, the dispersion relation will be divergent at the crossing point as well, and hence the 
perturbation will also not be stable. 

In summary, we have analyzed the most general case of a single scalar field described by a lagrangian in form of 
C = £(0, (9 M 0<9 M 0), and have studied different possibilities of w crossing the cosmological constant boundary. We 
showed that those cases can either cause the effective sound speed c 2 . to be negative, or lead to a divergent dispersion 
relation, which makes the system unstable. 

To conclude, we have proved that in FRW universe, it is impossible for a single perfect fluid or a single scalar field 
minimally coupled to Einstein Gravity to have its equation of state w crossing the cosmological constant boundary. 

To realize w across — 1, one must introduce extra degrees of freedom or introduce the nonminimal couplings or 
modify the Einstein gravity. In the recent years there have been a lot of activities in building models with w crossing 

-1 [aElzilil- 



The simplest Quintom model is to introduce two scalar fields with one being quintessence-like and the other 
phantom-like [16L l28l. l30f . However, this model suffers from the problem of quantum instability which is inherited 
from phantom [78l . l79l 180 . 181 1 . This issue could be solved in the effective description of the quintom model with the 
operator 0D0 involved. Consider a canonical scalar field with the lagrangian C — ^9^0(9^0— V((p). This model is well 
defined at quantum level, however, it does not give w crossing —1. As an effective theory as we know the lagrangian 
should include more operators. If these operators are functions of only the scalar field and its first derivative 
9^0(9^0, as we proved above the w still can not cross over -1. However as pointed out in Refs. 27, 2^], the existence of 
the operator 0D0 makes it possible for w to cross over the cosmological constant boundary. Furthermore, at quantum 
level as an perturbation theory this effective model is well defined. The connection of this type of Quintom-like theory 
to the string theory has been considered in [82| and 83 1. 
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